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Abstract 


We study two kinds of Dirac sea effects on the ^Sq pairing gap in nuclear 
matter based on the relativistic Hartree approximation to quantum hadro- 
dynamics and the Gor’kov formalism. We show that the vacuum fluctuation 
effect on the nucleon effective mass is more important than the direct coupling 
between the Fermi sea and the Dirac sea due to the pairing interaction. The 
effects of the high-momentnm cutoff are also discussed. 
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Superfluidity caused by the pairing correlation between two nucleons with the linear or 
angular momenta opposite to each other is a key ingredient to describe quantitatively the 
thermal evolution of neutron stars and the structure of hnite nuclei. As a way of description, 
relativistic models are attracting attention. The origin of relativistic nuclear models can be 
traced back to the work of Duerr |^. Since Chin and Walecka succeeded in reproducing 
the saturation property of symmetric nuclear matter within the mean-held theory (MFT) 
with the no-sea approximation |^, quantum hadrodynamics (QHD) has described the bulk 
properties not only of inhnite matter but of hnite spherical, deformed and rotating nuclei 
successfully |^,^. These successes indicate that the particle-hole (p-h) interaction in QHD 
is realistic. However, various observables of nuclear many-body systems are sensitive to the 
single-particle properties around the Fermi surface. The knowledge of the residual particle- 
particle (p-p) interaction is required to describe them. Since this is still less understood, 
non-relativistic interactions such as Gogny force are used in the pairing channel in practical 
“relativistic” MFT calculations. 

The hrst study of a relativistic p-p interaction for the pairing channel in the nuclear 
medium was done by Kucharek and Ring P . They adopted, as the particle-particle interac¬ 
tion (upp) in the gap equation, a one-boson-exchange (OBE) interaction with the ordinary 
relativistic MFT parameters which gave the saturation under the no-sea approximation. 
The resulting maximum gap was about three times larger than the accepted values in the 
non-relativistic calculations MM |IO|,pT|. Various modihcations to improve this result were 


proposed. They can be classified into two groups: One is to adopt the Upp which is consis¬ 


tent with the p-h channel 


and the other is to adopt effective ones which are not 
Among the former, we examined the 


15,3,16 


explicitly consistent with the p-h channel 

p-h polarization in Upp which reduced the pairing gap in the non-relativistic models |1l7|,p!8 


But the result was negative; this suggests that the nucleon-antinucleon (N-N) polarization 


should be taken into account simultaneously |T^. Before doing this, the antinucleon degrees 


of freedom have to be taken into account in the OBE step as did by Guimaraes et al ag¬ 
in this paper, we discuss a more important effect of the Dirac sea than that they discussed. 

Our formulation consists of three steps. In the hrst step, the equations of motion of 
the normal and the anomalous Green’s functions, G and F, are derived. We start from the 
ordinary a-u model Lagrangian density. 


C = — M)ijj 


+ Qa'ilJcrilj 






dyUJ^. 


( 1 ) 


The Hamiltonian density is derived from C and the chemical potential p is introduced in 
the Hamiltonian, H' = JiTi. — The equations of motion for 

Gapix - x') = -f( 0 |T?/)„(a:)^/ 3 (a;')| 0 ), 

Fa,p{x - x') = -z(0|T?/>o(a;)?/5/3(a;')|0), (2) 

are derived from = [ijj, H']. Here we note that the Green’s functions have to be dehned 
by the superhuid ground state |0) in order to introduce the antinucleon in the next step. 
Aside from this, this hrst step is essentially the same as the formulation of Ref. [Q. After 
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evaluating the commutator with the interaction term in H', the meson helds are eliminated 
by using the inverse of the Klein-Gordon equations. Here the gauge term in the propagator 
of the LJ meson can be discarded since to couples to the baryon current which is conserved 
in average in the present superfluid case. The result is 


M + - x') = - x') + y)'tp 5 {t, y)^p^{x)^jJp{x')\0), 

{-p - M + 'y'^fi)^aF^i 3 {x - x) = 

'^a'y,&e QaFcr^a'y^Se T 


with X = {t, x) and p = i(3{dt — a ■ V). Here we introduced an instantaneous approximation 
since it was reported in the preceeding works [^,^] that the retardation effects are small. 
The time-ordered products in the right-hand side are decomposed a la Gor’kov |T^ by 
extending Wick’s theorem, 


y)Mt, y)^p^ix)'^fs{x')\ 0 ) 

= -peS'iOhs'S^G'ipix - x') + p.^5/(x - y)ylsiGe 0 {,y - x') - «:^,(x - y)iFs 0 {y - x'), 
(0|TV',(t, y)pJ5{t, y)i)^{x)Pl) 0 {x')\Q) 

= -peS'i^hlsiF^yix - x') + p,y(y - x)7°yzF5^(y - x') 

+«vy(y - ^)ll5ly^'^G^0{y - x), (4) 

where the normal and the anomalous densities are 


Pa/3(x-y) = (0|?/^J(f,y)^/>„(x)|0), 

«:a/ 3 (x-y) = (0|V’/3(t,y)^a(a:)|0). (5) 


The selfenergy and the pairing held are dehned by these densities as 


Eo^(x - y) = ( 14 ^ 7 ePe 5 '( 0 ) - 14 , 7 ^p, 5 /(x - y))7°/5, 

A„^(x - y) = 7L'7?y''4'^.7'^«<5^(x “ y)- (6) 

The Fourier transform of the equations of motion for g and / dehned by = 5 ' 7 / 3 ' 7 y^ and 
F ^0 = /y/3'7y77y/3 is given by 


( h /r)a7 ^a'f ^ f 9'y0 A _ f ^a0 ^ /y\ 

where h = a • k -|- f3{M -|- S(A;)). 

In the second step, we derive the equations for the Bogoliubov transformation amplitudes 
by expressing g and / in Eq.(|^) in terms of them. We introduce the Dirac held in the 
Schrodinger picture. 


^(x) 


i: (aX(A)e““ 


+ 6iF(A)e-‘*^"), 


( 8 ) 
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to manipulate the t-dependence. Here the normalization of U and V is chosen conforming 
to Ref. 


23 . The Bogoliubov amplitudes are dehned by 


Ax = ( 6 |oi>)j,| 0 >, Bx = ( 6 | 6 LaI)a| 0 ). 

Ca = (6|oLa<;1|0>. £>a = (0|6ai71|0>, (9) 

as the overlaps between the quasiparticle J]x> 




( 10 ) 


and the nucleon or the antinucleon. These give the explicit expressions of g and /, 

g-,s{cj,k) ^ —l—Y,{\Ax\^U-,{\)U;{X) + |BaPU(-A)V';(-A) 

^ "r 5 

+AxB-xU-,(\)V^{-\) + BxA-xV-,(-X)U;m 
+ . i:{|g-ArcfA(A)i 7 ^(A)+iD_Aru(-A)i/;(-A) 

^ “T C}^ It g 

+D_xClxUx{X)V;(-X) + C-xD'_xVx(-A)U;{\)}, 
f., 0 {u,k) - ,_y ^■ Y.{c>.AXu;(-x)u;{\) + DxB‘xV;{x)v;{-x) 

+CABiy (-A)VJ(-A) + DxAlV.J(\)U‘0{X)} 

+--T-rzT,{C->.A‘_xU;{-X)U;(X) + D,xB‘_xV;{X)V^(-X) 

+B_aV Ay (-A)Vg(-A) + C-xB'_xV.:(X)U^{X)}, ( 11 ) 


the former describes the normal propagation and the AN = 0 pairing, and the latter de¬ 
scribes the ordinary AN = 2 pairing. This AN = 0 pairing appears because Cooper pairs 
can be formed regardless of the sign of the single-nucleon energy. Substituting Eq. ([TI|) into 
Eq.(^ and dehning the matrix elements of the pairing held A as 

A(A) = 

A(A) = e-'“-(-"Vi(-A)A„^T,;3R^(-A), 

5(A) = e-*“-(-")Ht(A)A„.,T,^l/^(-A), (12) 


for the Fermi sea pairing, the Dirac sea pairing, and the AN = 0 pairing, respectively, we 
obtain the equation of the Bogoliubov amplitudes. 


— Ek + g, 

0 

-A(A) 

-5(A) \ 


/Ax\ 

0 

u! -|- Ek -l- g 

-5(A) 

-A(A) 


Bx 

-A(A) 

-5(A) 

UJ + Ek- g 

0 


Cx 

-5(A) 

-A(A) 

0 uo 

Ek J 


\DxJ 


Here we used 

ha^U-y{X) = EkUa{X), ha-yV^{ — X) = —Efcl4( —A), 
Ek = Vk2 + M*2, 

TU{X) = e*“+(^)l7*(-A), TV{-X) = e*“-(-^V*(A), 
Th* = hT, TA* = AT, T = 


(13) 


( 14 ) 
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Note that we adopted notations such that Eq.(|I^ took the same form as that of Ref. 
and chose a phase convention such that all the matrix elements were real. Among the 
eigenvalues, 


a;" = {El + + + A^) + ^2) _ - A2 + + 4d^{4El + (A + A)2) (15) 

corresponds to the Fermi sea pairing in the decoupling (h —>■ 0) limit. 

In the third step, we express k, and subsequently the matrix elements A(A), A(A), and 
A (A), in terms of Ax - Dx- The Fourier transform of k in Fq.(D is given by 

Ks,{uj, k) = 27rA(a; + 4) E {A-xC.xU,{-X)U,{\) + B_xD_xV,{X)Vs{-X) 

S 

+A.xD.xU,{-X)Vs{-X) + B_xC-xV.{X)Us{X)}. (16) 


Substituting this into the Fourier transform of the second equation of (^ gives 

A(A) = 4 

[Tr{(^ + M*){^ + AxCx + Tr{(^ + M*)(^ - BxDx] 

, 9l 


r dXp 1 / 


9l 

/ { 2 TifAEkEy 

k — p 

p + 


|k — p|2 + 

[Tr{(^ + M*)7 Ai^ + ^*)7/.}e*“+^^^AAC'A+Tr{(^ + M*) 7 A^ “ BxDx]), 

A(A) = 

[Tr{(^ - M*){i) + M*)}e*"+(^)AACA + Tr{(^ - M*){f> - BxDx] 

, al 


f d?p 1 ( 


9 I 

1 {2n)^4EkEpy 

■|k-p| 

p + ml 


|k — pP + 

[Tr{(| - M*)7 Ai^ + M*) 7 M}e*“+^"^AlAC'A+Tr{(| - M*)7A^ “ BxDx]), 

6 {X) = -- 


CO 

— 

9 I 

1 (27r)3 4FfcFE 

|k-p| 

p + ml 


[Tr{75(^ - M*)7V(;^ - M*)7°}e“+(^)RAC'A 
-Tr{7''(^ - + M*)7°}e'“-^"^^AADA] 


|k — pp + 

[Tr{7A^ - M*)7°7 m7'(|^ “ BxCx 

-Tr{7''(^ - MA7°7 m 7''(^ + M*)7°7M}e“-^""^7lADA]), 


(17) 


where both 7 ^ in each line are covariant due to the time reversal. Here we used A_xD_x = 
—BxCx, and so on, derived from the dehnition of k, the s-independence of e®“+*^^MAC'A, and 
so on, and 

j:u(\)m) = Nr(-A)i/(-A) = - An, 
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( 18 ) 


i:C/(A)i/(-A) = - M*)7», 5 ;: V(-X)U{X) = + M*)7». 

k = {Ek, -k). 


Equations (|^) and (p!7D form a system of selfconsistent equations. The nucleon effective 
mass M* is determined by the selfconsistent condition for the scalar part of the selfenergy, 
the hrst equation of (0), 


M* = M 


9o 


9 a 7 


1-Ac 


M* 


cydp 


m: 


27r^ JO a/p^ + M*'^ ^ 


mi TT^ 


1 M* "a 11 

{M*3 In (—) - M\M* - M) - -M(M* - Mf - —{M^ - Mf}, 


(19) 


within the relativistic Hartree approximation (RHA) in which the divergence due to the 
Dirac sea is renormarized using the counter terms 


/let (y\(7 




( 20 ) 


Hereafter we suppress the argument s, and therefore the direction of k which dehnes the 
direction of s, of the gaps and the Bogoliubov amplitudes because they affect only the 
overall sign. The expression of the vacuum contribution in Eq.(]^) is taken from the non¬ 
superfluid case 1^. We will discuss this later. Then the actual task is to solve the coupled 
equations (11^), (pTl), and (^). If we neglect the vacuum fluctuation contribution in Eq.(^), 
the system of equations corresponds to that of Ref. [|^. Here we note that we adopt the 
Hartree approximation with or without the Dirac sea contribution as in Ref. for the 
no-sea case. 

Now we proceed to the numerical results. Parameters used are M = 939 MeV, = 550 
MeV, = 783 MeV, g'^ = 62.89 for the RHA or 91.64 for the MET, g'^ = 79.78 for the RHA 
or 136.2 for the MET [^, and 7 = 4 (symmetric nuclear matter). First we look into the 
relative magnitudes of the gaps A{k), A{k), and 6 {k). Since Aj. and are dominant among 
those associated with the eigenvalue Eq.(^), A{k) and A{k) are much larger than 5{k). 
From the dehnition of n in Eq.(|^), two kinds of terms, one is antisymmetric in spin space 
and even with respect to the inversion of the momentum, and the other is symmetric and 
odd, are possible in K 5 e(a;,k). If the hnite-range effects are neglected, only six terms of the 
former type are possible as discussed in Ref. |^. They correspond to scalar, pseudoscalar. 


and 4-vector terms in kT. Among them, the scalar and the time component of the vector are 
dominant and consequently A{k) and A{k) are determined by these two terms as pointed 
out in Ref. |T^. Basically this applies also to the calculations of Ref. |T^ and ours in 


which the terms of the latter type are also included. In contrast, 5{k) which measures the 
AV = 0 pairing is given by the pseudoscalar and pseudovector terms in kT. Since they 
contain B^Ck and A^Dk, their typical values are the order of 10“® MeV or less both in the 
RHA and in the MET. Therefore, as for the magnitude of the AN = 0 coupling effect on 
A{k), our calculation does not agree with that of Ref. [^|. The origin of this disagreement 


is that their Eq.(54) does not have this structure that the different type of products of 
the Bogoliubov amplitudes appear in 6 {k). The present result indicates that the Fermi sea 
pairing and the Dirac sea pairing decouple in accuracy of 10“^. Accordingly, we use the 
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expressions for 6 {k) =0 hereafter. Since the eqnations for Ak and Ck deconpled from Bk 
and Dk are eqnivalent to a gap eqnation, 




p^dp, 


( 21 ) 


where Vpp{k,p) is an antisymmetrized matrix element of the adopted p-p interaction with 
an instantaneons approximation and an integration with respect to the angle between k and 
p to project ont the 5'-wave component, the nnmerical task is greatly simpihed to solving 
the conpled eqnations (|T^ and (^). 

Next we compare the resnlts of the RHA and the MFT obtained by adopting Ac = 15 
fm“^ which is large enongh for the nnmerical integrations to converge (see Fig.^). These are 
presented in Fig.Q. This shows that the Dirac sea contribntes to redncing the gap at low 
density while to enhancing it at high density. This can be nnderstood as follows: Althongh 
both Qu and are reduced in the RHA in comparison with the MFT, altogether they act to 
reduce the repulsion as shown in Fig.Q. Since both the low-momentum attraction and the 
high-momentum repulsion give positive contributions to A{kp) as discussed in Refs. [P!5| , P!B| , 
the above-mentioned reduction of the repulsion leads to reduction of A(/cf) in the RHA. 
Besides, the low-momentum attraction decreases steeply as density increases especially in 
the MFT (not shown). This leads to steeper reduction of A(fcF) at high density in the MFT 
than in the RHA. Here we note that this RHA calculation was done by assuming the same 
vacuum fluctuation contribution as the non-superfluid case as mentioned above. We believe 
that this approximation is practical because the effect of the pairing on M* is negligible 
numerically except at very low density. This indicates that the bulk property, M*, affects 
the Fermi-surface property, superfluidity, whereas the opposite is not true. 

Since QHD is an effective theory of hadrons, form factors or cutoffs related to the spatial 
size of hadrons might be necessary . On the other hand, an important feature of the gap 


equation is that it has such a form that the short range correlation is involved |1^ , ^ , |9[| . This 
leads to a similar oscillatory behavior of A{k) to that of Vpp{k, kp) as functions of k 
Matera et al. found that their formulation for the pairing gap gave a cutoff 


15,16 


13 . This is 


interesting in respect that both which are originated from the pairing correlation act to evade 
the repulsion at high momentum. They reported that the values they obtained were 1.7 - 
1.9 fm“^. These values mean that the repulsive part was completely cut, see Fig.g Since 
our formulation does not have any means to choose a cutoff as those of Refs. p|,|^ , here we 
show the dependence of A(A:f) at k-p = 0.9 fm where it becomes maximum, on the cutoff 
momentum Ac as a free parameter. The result is presented in Fig.^ This shows that the 
MFT result decreases steeply as the cutoff decreases. This can be understood as follows: 
The difference between the RHA and the MFT in the low-density case in Fig.|^ mainly 
comes from the difference in the magnitude of the high-momentum repulsion. In other 
words, the contribution of the high-momentum repulsion is more important in the MFT 
case. Consequently, as the cutoff momentum decreases, the A(A:f) of the MFT calculation 
decreases more. The plateaus around Ac = 2 - 3 fm“^ in the RHA and Ac = 1 - 2 fm“^ 


m 


the MFT are due to the sign change of Vpp around there (see Fig.^. A comparison with 


the Bonn-B potential which reproduced the maximum pairing gap accepted in the non- 
relativistic studies ||I5|,|3| is also shown in Fig.^. This indicates that the high-momentum 
repulsion of the MFT without a cutoff is too strong and the low-momentum repulsion is 


7 











































too weak both in the RHA and in the MFT to give reasonable pairing gaps. A possible 
improvement in the p-p channel would be to take into account the nucleon-antinucleon 
(N-N) polarization. 

To summarize, we have studied two kinds of Dirac sea effects on the pairing gap in nuclear 
matter. One is the vacuum fluctuation effect on the nucleon effective mass and the other is 
the direct coupling between the Fermi sea pairing and the Dirac sea pairing. The former is 
a bulk effect while the latter affects only the Fermi surface. Our calculation indicates that 
the former is more important. The dependence of A(/cf) on the high-momentum cutoff was 
also discussed and the MFT result has been shown to depend more strongly on it than the 
RHA one. 

Discussions with Prof. M. Nakano and Prof. R. Tamagaki are acknowledged. The author 
also thanks Mr. T. Tanigawa for calculating the Bonn potential. Numerical calculations were 
done using the computer system of the Information Processing Center, Fukuoka University 
of Education. 
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FIGURES 


Pairing gap at the Fermi surface 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

kp (fm-i) 

FIG. 1. Pairing gap in symmetric nuclear matter at the Fermi surface as functions of the Fermi 
momentum. Solid and dashed lines indicate the results obtained by taking and not taking into 
account the vacuum fluctuation contribution, respectively. 


Particle-particle interaction 
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FIG. 2. Matrix element Vpp{k,k-p) as functions of the momentum k, with a Fermi momentum 
kp = 0.9 fm“^. Solid and dashed lines indicate the results obtained by taking and not taking 
into account the vacuum fluctuation contribution, respectively. The Bonn-B potential for the 
channel is also shown by the dotted line. 
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Cutoff dependence 
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FIG. 3. Pairing gap at the Fermi surface, /cp = 0.9 fm“^, as functions of the cutoff parameter 
in the numerical integrations. Solid and dashed lines indicate the results obtained by taking and 
not taking into account the vacuum fluctuation contribution, respectively. 
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